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TWO-DIMENSIONAL POTENTIAL FLOW PAST 
AN ORDINARY THICK WING- PROFILE* 
By F . Keune 

} 

SUMMARY * 

This report deals with the development of a method 
which gives a lucid and convenient solution of the flow 
conditions in the vicinity of a common, thick airfoil sec- 
tion wherein the thickness of the profile is taken into 
account. The method consists in making the airfoil the 
streamline in a parallel flow hy disposing on its mean 
line certain source and vortex distributions the fields 
of which are superposed on the parallel flow. These dis- 
tributions of singularities are secured for the general- 
ized Ka'rman-Treff tz profile by means of conformal trans- 
formation from the flow about a circle. Five different 
distribution functions are afforded for the density of 
superposition, which combine in a specified manner to the 
necessary distributions of sing\ilarity and represent a 
generalized Kar man— Tref f t z profile in parallel flow. For 
these profiles the speed for each of the five distribu- 
tions is then computed independently of the angle of at- 
tack . 

This substitution of a profile theoretically secured 
by singularities is then extended to include more general 
airfoil shapes, since a simple method permits a close ap- 
proximation of any profile by a generalized Kar man-Tr ef f t z 
profile with known substitution of singularities. This 
approximation then affords the solution of the velocity 
field in the neighborhood of the general thick airfoil by 
a convenient method, For each of the five functions of the 
superposition density a speed diagram is plotted; the 
curves of constant speed parallel to the coordinate axes 
are c omb ined„ t o an orthogonal system. The readily dis- 
cernible effects of the individual singularity distribu- 
tions on the speed at any point of the plane of the section 

*"Die ebene Pot ent ials tr omung um allgemeine dicke Tragflug- 
elprofile." Jahrbuch 1938 der deutschen Luf tf ahrtf or schung , 
pp. I 3 — I 26. 
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are in a simple manner correlated to the speeds. For com- 
parison the speeds can he computed by the same method con- 
formally to Birnbaum's approximation. It is found that 
the omission of the profile thickness is permissible only 
at great distances , that is , approximately of the order of 
magnitude of a wing chord. 

A method is also given which affords a check on the 
agreement of the profile contour with the closed stream- 
line defined by the singularity distributions. A further 
method enables the prediction of the profile contour for 
any selected singularity distribution. 

After completion of this article, Z.f.a.M.M., vol. 18, 
1938, p. 107, carried a report by P. Weinig, entitled: 
"The Flow Conditions in the Pield of Thin, Slightly Cam- 
bered Airfoil Sections," which gives the velocity distri- 
butions for thin profiles in accord vrith the writer's 
calculations (figs. 26, 27, and 28;. V/einig also consid- 
ered the effects of complicated camber shapes. But the 
essential of the present article is the determination of 
the effect of profile thickness, which was not considered 
in Weinig's report; although it is very pronounced, accord- 
ing to present study. In general, it' should be stronger 
than the effects of the mean— line forms exceeding Birnbaum's 
f or aulas . 



IHIEODUOTIOH 



A solution of the two-dimensional flow about thin 
airfoil sections is afforded from either Birnbaum's method 
(reference 1) or G-lauert's (reference 2) by the superpo- 
sition of vortices on the mean line of the profile. But 
this theory fails to allow for the profile thickn_e_ss. 
While the two-dimensional flow about thick profiles can 
bff solved by means of conformal z ransf or mat; i on , it in- 
volves considerable paper work, except for particular 
families, such as the .Joukowski, Kar man— Tref f t z , and gen- 
eralized Xarma'n— Tref f tz airfoils (refere nee 3 , . In addi — 
tion to that, these arguments do not lend themselves to 
the solution of the flow in vicinity of a wing of finite 
span, because conformal transformation i_s a pplicab le only 
to two-dimensional flow and the application of the dataT 
derived from it to three— dimensional- processes is limited 
to particularly favorable cases. On the other hand, these 
interference fields in three— dimens i onal flows are fre— 



NAC.A Technical Memorandum Ho. 1023 



3 



quently important, especially for "biplane calculations,* 
Kence if .the. thickness effect of a profile in all such 
calculations is to "be accounted for, the conformal trans— ^ 
formation must be abandoned in favor of the singularity 
method as for the thin profile. 

The method of singularities consists in visualizing 
within the body in a stream ano.ther stream that forms the 
analytical continuation of the stream outside the body. 
The total flow has a closed streamline which coincides 
with the body contour. ?or the flow within the analytical 
continuation yie.lds singularities, namely, sources and 
vortices. The total flow can be interpreted as the field 
of these singularities. Thus the plane flow about a cyl- 
inder is represented by a doublet superposition on a par- 
allel flow (that is, the flow from a source to a sink lo- 
cated infinitely close to one another). The cylinder 
radius is defined by the ratio of doublet intensity to 
speed at infinity. If the cylinder has lift, a vortex 
flow past the cylinder axis must be added. In general, it 
requires several, usually continuously distributed singu- 
larities rather than singularities in a single point as in 
this instance. On surfaces of the type of airfoil sections 
the singularities are usually located along a line, that 
is, on the mean line. The displacement of the streamlines 
from the mean line corresponding to the finite thickness 
is e.ffected by a source distribution; whereas a vortex dis- 
tribution produces the lift-inducing def lect ion . of flow. 
Such representation of flow past a profile is always pos- 
sible if the flow in the outer space can be analytically 
continued on the inside (theory of analytical functions) 
so' that the singular points of the function are disposed 
on a line wholly within the profile. This is necessary 
because it is conceivable that isolated singular points 
might yet be present; although they are likely tc occur 
only along the chord by very radical thickness changes, 
such as are not apt to be encountered on the conventional 
airfoil shapes. 

The principal task of the present report is to evolve 
a simple method by means of which the intensity and the 
position of the singularities can be so defined that the 
closed streamline resulting from the super pos it ios of the 
singularity' flow on the parallel' "f low is" coincident with 
the given profile contour. 



*The solution of this' problem forms the subject of a spe- 
cial report. 
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If a symmetrical profile is chosen, the singularities 
will, for reasons of symmetry, he located on the mean line. 
Then the singularities could he determined in a similar 
manner as known from the calculation of the flow past air- 
ship hulls, hy assuming the source distributions partly 
constant as proposed by Von Karman (reference 4) and so 
defining the unknown intensities that the hody contour he- 
comes stream surface. But then the determination of the 
starting point of the superposition involves unusual dif- 
ficulties. If it is wrongly chosen, the density of super- 
position at the start fluctuates "between positive and 
negative values, which was the very reason for the ulti- 
mate change to surface superposition on airship hulls, 
especially in the case of transverse flow so as to avoid 
these difficulties (reference 5). For the two— dimens i onal 
problem in question a different line of attack is pursued 
from the very beginning. Since the flow past airfoils of 
a certain family is readily afforded from the flow past 
a circle by means of conformal transformation, the singu- 
larity substitute is forthwith computable for such pro- 
files after inclusion of the flow xvithin the conformal 
circle. Certain laws and typical distribution functions 
are acquired here which are faithfully transformed into 
general profile forms. 

If the singularity distributions are given, the next 
problem consists in computing the flows of the separately 
found singularity distributions which are then easily com- 
bined into the desired flow of" a prescribed profile. 

With the computed singularity distributions new pro- 
file forms with desired form characteristics can also be 
constructed without having to resort to conformal trans- 
f ormat i on . 
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y^ ordinate of the streaaline 

a * 

y x - ordinate of the mean M-ne at— point- x - + — 
y s ordinate of the mean line at point x = — — ' 

Di thickness of profile at point x * + - 

2 

D 2 thickness of profile at point x = — — ' 

2 

D max mft 3cimum thickness of profile 
t profile chord 

a distance OA 2 , criterion in plane of circle 

a 1 half distance of "branching points in plane of section 

d measure for profile thickness 

tan 0 measure for profile career 

H measure for reversed curvature of profile 

6=ktt trailing edge angle of profile in radians 

$(cx), $(z) complex potential function in plane of circle 

and section, respectively 

^(tf). (z ) stream function in plane of circle and sect i on , 

respect ively 

Vco e~ ia ' air flow in plane of circle 

Vco e air flow in plane of section 

V x s=Vco cos a component of flow in x direction - 

V y=Vco s * n component of flow in y direction 

v n normal component of velocity at the slit in plane £ 
and z, respectively 

v^ tangential component of velocity at the slit in 
plane £ and z, respectively 

u speed along x axis in plane of section 

v speed along y axis in plane of section 
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w absolute amount of speed in plane of section 

* angle of speed toward positive x axis of plane of 
section 

/ f > 

^v^' ^vy^T ^ strength of source per unit distance 

Y v ( ! strength, of vortex per unit distance 

/ . \ Va ' 

f v I £7 j "base function pf.the singularity distributions 
V / with the ordinal number u 

/ x y N 

Qui —J -7 ratio of integral of flow function of source 

\ a a V , * s 

distribution f ' _ ; to undisturbed 

V a y _ 

air flow at point --; 

a ; a 1 

( x y \ 

w u 1 ~ ■ ratio of the integral cf the flow function 

f 1 \ 

of the vortex distribution f v ( -2- 

\ a 1 / 

I. THEOHY 
1. Solution of Singularities 
a) J^ukowski Profile 



The Joukowski profile K s 1 (fig. 1) in the £ plane 
in the reflection of a circle If 3 in the a plane, 
transformed by the conf ormal function 

I = o + £ (1) 

frome plane a in plane £. This circle K 3 is given 

by three quantities, namely, distance AjAg = 2a, \vhich 

defines the profile chord, distance 00 z = a tan 
which defines the profile camber and distance 0 X 0 B = d, 
which defines the profile thickness. Prom these quanti- 
ties the radius of circle K 3 follows at 

R = — ^— + d (2) 
cos P 

its midpoint 0 e has the coordinates 



a Q = - d cos 0 + i (a tan 0 + d sin 0) 



(3) 
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Circle K,, with center 0 lt and radius and a 

cos P 

radius vector amounting to 1 ' 

r ~ a .1 + tan p sin tp + i tan 3 'P sin 2 "9^ (4) 

(fig. 1) is transformed by the c onf ormal .f unct ion (-1) in 
a double curve K : ' , considered herein as the mean line 
of the profile. Its coordinates are: 



£ £ /la 2 N 

— = — «s cos cp ! 1 + ~ tan"* p sin" cp ] 

2a a' V 2 J 



s 
1 
1 



(5) 



-- = --« tan P sin 8 9 - W riO - cC<p\ 
2a r ' 



?»r 



-2- * tan 0 (l - -L-^. (6) 

? . 

The slope of the mean line at any point — - is accord- 
ingly a 

dn 

— = tan v « - 2 tan p cos cp (7) 

The inside of this circle K 1 then forms the second posi- 
tion of the I plane, so that the mean line is at 
the same time the slit in the original position for tran- 
sition into the second. Points A, 1 . and A 2 ' are the 
"branch points of the transformation. The_distance of these 
points from the zero point is U T A"7 T = U'Ag"'' = 2a = a'. 



The flnw past profile Kg 1 is obtained by transform- 
ing with function (l) the flow past circle K 2 from plane 
a into .plane £ . This f low in pla ne o , wherein circle Kg 
is in a parallel flow at angle a 1 .and speed Vo, at in- 
finity, is given by the complex flow.' potential 

. B s la' , r 
4(a) = V te j (a-a 0 )e ia + + i — In — 2. ( 8 ) 

If, to insure smooth flow— off at the trailing edge, a stag- 
nation point is to be located in point A lP then it must 
be that 
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r = 4tt 7co & sin (a • + P) (9) 

Now when the flow conditions within the^'pr of ile "are 
considered as continuation of the looked— for flow in the 
outer space, it is hardly to "be expected that the same 
speeds prevail on the upper and lower edge of the assumed— 
ly double slit %-x', because these edges form the transi- 
tion to the. second position. This is .readily apparent 
upon examination of the speeds at circle Ej of plane a 
the upper arc AjBjAg of which "becomes the upper edge 
and the lower arc A 2 B 2 A i becomes the lower "border of the 
slit. The flow past circle K 3 affords the speed at 
circle Z 1 ; whereas conformal function (l) yields the 
speed at the slit which can he divided in its normal and 
tangential components. To get rid of the conformal trans- 
formation, the plane of the profile, plane £ , must he 
interpreted as plane closed in itself, wherein the slit 
loses its significance and becomes the mean line. Then, 
however, the difference of the normal speeds v Q on the 
upper and lower border of the slit must be interpreted 
as source superposition (see fig. 2) and that of v^. as 
vortex superposition. The s ource' strength q. per unit 
length at a point ^ of the mean line is accordingly 
given by 

4( k) = <v n U - (v Q ) t (10) 
n ob ^1 a u i-j 

and the vortex strength V per unit length by 

With these data the singularity distributions for a 
Joukowski profile can be computed. 

For future purposes the flow velocity in plane 

cr (fig. l) is divided in its components parallel to the 
axes of the coordinates 



v w = V oo c os a 1 
V_ = T w sin a 1 



(12) 



i 



and each of the flows created by the components is analyzed 
separately. The complex flow potentia-1 (8) divided con- 
formally to (12) affords the speeds v CT in the a plane 

and. the speeds v,£ = u — iy through the known relation 
v^ = v^ ™. The speed in the £ plane produced by the 
0) flow reads 
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d 4^ 'd $ (ct)^ do- 



d a 



4t 



a 
a 



R 3 fa 



a 2 Lj 



-- + 3^- — 

_ 2 S 3 

a a a cy _ 



R r a 0 a 

+ 21- sinp 1 1 + — - + 



cr 0 2 a 3 



a 



a 2 J 



(13) 



a 
a 



a 
a 



and that due to the t flow is: 



d $ T d $ (a) T da 



a R 2 



7, 



d a 



a o. 
-+2 — 



at 



a 0 2 a 3 
+3 _£ 

a a 



R 



- 2 - cosp Jl + 
a La 



a 0 a 



- + — 

_ _ 2 _ 2 

cr a a 



^1 

_ s 



- (14) 



a 

a 



a 
a 



on the assumption that the distance [ff 0 | of the center 
0 2 of the circle from the zero point is short compared 
to a and hence small compared with \a\ = r. This is 
always the case when restricted to the conventional profil 
where d/a is always <K 1 and tan p 1 ; hence, es- 
pecially thick strut sections are discounted. 

Entering the values (2), (3), and (4) in equations 
(13) and (14), where a = r e ic P indicates the coordi- 
nates »f circle K lt gives the velocities on the mean 
line, the real and imaginary part representing the com- 
ponents u and (— v) al^ng the coordinate axes. From 
these, with consideration to equation (7), the tangential 
and normal components follow at 



' t = u(l — 2 tan 3 p cos s cp) -r 2 v tan P cos cp "j 
r n = 2 u tan p cos cp + v(l - 2 tan 2 P cos 3 cp) J 



r 



(15) 



(see fig. 3) . 



After expansion of the speeds (equations (13) and 
(14)) in powers of d/a and tan p — the terms of higher 
than the second order are considered negligible — there 
is obtained with equation (15) and consideration of (10) 
and (11) a close approximation for the singularity dis- 
tributions of the Joukowskl profile. 
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For the flow along the u> axis the source distribu- 
tion reads* 

q.(ep) = - 2 V w — s -[cos 9(1 - cos 9) + sin 2 cp] I 
sin (p (a L T J I 

^2 

+ 2- — [cos cp(l - cos cp) + sin 2 cp - 3 sin 2 cp cos cp] + . . . > (16) 
a j 

and for the vertex distribution 



y (cp) = 2 T a> <^2 tan P sin (p 

and for the source distribution along the t axis 



d 1 — cos cp 

- - tan 3 | — - 4 sin cp 

a L sin cp 

1 

+ 6 sin cp cos cp j + . . X (17) 



q. (cp) = 



,, d n cos cp(l - cos cp) 

-6 V T - tan P j - — + sin cp 

' a (_ s in cp 



- 2 sin 9 cos cpj + ... (18) 

and the vortex distribution 

/ \ ~ „ fl — coscp _ d,l — coscp ~|| 
Y (cp) = 2 V T J +2 - i - sin cp i I 

\_ s in 9 a Lsincp Jl 

+ 3~i ±_r_£££_c£ 2 sin 9 + 2 sincp cos q, I 
a l s in cp J 

— — tan P i — ° ° S ~ — 4 sin 9+6 sin 9 cos 9 1 

a L sin 9 ' j 

S(l"*"COSCPl "** 

- tan P I sin 9 + 2 sin 9 cos 9 + . . .± (19) 

L s in 9 2 j J 

Between the angles 9 of plane a and the running 
coordinates of plane £ the simple relation 



*The heavy separation lines in the formulas indicate the 
places vhere the equations should be broken off when terms 
of the second degree are discounted. 
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cos 



tp K - 

a' 



(20) 



exists according to equation (5); besides , the flow at 
infinity remains unchanged in quantity and direction in 
the transformation, so that = V x , V T = V , and 

a 1 = a . 

A study of the computed singularity distributions 
discloses a consistent recurrence of the individual ana- 
lytical expressions for the superposition density from 
which the entire distribution can be obtained. They are 
the four analytical expressions 



f i(£/a') = — [cos cp(l - cos cp) + sin 2 cp] 

sin cp 



f a(|/a') 
f 3 (£/a«) 



+ /h~±i±l (1 + 3 e/al) 

v 1 + |/a' 

AT"T7a T 



= + 



a 

1 — cos cp 

sin cp J 1 + | / a 

sin cp = + y 1 -{^ j 



> 



(21) 



f 4 ( |/a 1 ) = sin cp cos cp = — /l 

a 1 v 



which, for practical reasons, are treated individually. 
The expression for the source distribution f 1 {|/a 1 ) is 
given separately; although it is composed of f 2 and f 3 



f !(|/a') = 2f 3 (|/a') - f a (|/a •) 



(22) 



because of the necessary secondary condition 

j f ft/a')d| » 0 (23) 



+ a • 

_i_ r 

a~ 



—a 
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which must accompany all source dis tr ibut i oils . Through it 
the contour stipulated "by the source distribution 

/ t > 

f \^~~ J "becomes a closed curve at infinity. Functions 



t,(±) and tl±) 



\ -a- ■ ana i 4 i i of which the source distributions 

V a ' J \ a ' J 

(16) and (18) are composed, each satisfies this condition 
(23) "by itself. 

In the distributions f 2 (~ ; and fs/'--^} there 



can be recognized the well— known Birnbaum formulas (refer- 
ence l) for the vortex distributions which transform a 
flat or curved plate into a streamline. It affords, ac- 
cording to (19), the vortex distribution for the plate which 
is always occurring 

Va» / 7 V i + |/ a . * V a J 

and, according to (17), the vortex distribution which indi- 
cates the effect of the curvature and produces a rotation 
of theizero lift direction (a 0 = — P ) : 



\ 



Ys^^-y = 4 V x tan & Jl - (|/a') 2 = 4 T x tan 3 f 3 

(25) 

These known Birnbaum distributions are supplemented 
further by the effects due to profile thickness, that is, 
the source distribution (equation (16)): 

* 1 V a ' / x a7 1+ ^/ a ' V a'/' Xa Va » J 

(26) 



and the vortex distribution (equation (19)) 
Y23 



W 7 a v ' 1 + e/a , a , 7 a u * 3 Va V % a 



e/i 

(27) 

The quadratic terms are discounted for the present. They 
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merely afford in closest proximity of the profile contour 
a portion to the fow and are needed only when checking the 
-coinc idence - of --t-he— c-l-os-ed «s -or eamline resulting from the 
singularities with the profile contour (cf. appendix). 
Otherwise they are important only for error appraisal. 

, All singularity distributions extend from | = +a ' 
to | = —a 1 , that is, from the trailing edge of the pro- 
file, the rear branch point A a « to the forward branch 
point A 2 ', located slightly "behind the nose. The dis- 
tributions f i ( --^ and f s ( have in the forward 

<a 1 /' \ a 1 / 

branching point |= —a 1 an infinity point, that is, they 

N 1 

approach infinity as • ±— > s 



b) The Karman-Treff tz Profile 

Whereas the Joukowski profiles terminate in an infi- 
nitely thin trailing edge, the Karman-Treff tz profiles 
form a finite angle 8 on the trailing edge and so ap- 
proach the commonly employed profiles mush closer. Now 
the question arises as to the changes in the established , 
distributions when a Karman— Tr ef f t.z profile is to be stream- 
line in a parallel flow. The conformal function 



l_ = _ak 
1 + ak 




(28) 



transforms the circle K x (fig. 4) in a crescent K 2 ' 
with angle 6 , when 

2tt - 6 . 
k = = 3 — K (29) 

TT 

and 

K = (30) 

For .6 = 0, circle K a becomes the infinitely thin mean 
line, of the Joukowski profile, since (28) affords the 
conformal function (l). The Karman-Treff tz profile will ^ 
be, on the whole, thicker by the thickness of the crescent 
than the Joukowski profile. In view, therefore, of the 
singularity distributions a greater disp_lacement flow may 
be_ anjy^cipa^ed , that is, a new__ source distribution; where- 
as the^ vortex distri buti ons change "very lit.kLe." Our" 
^t^lysls -will proceed theref orVfron '» symmetrical 
Aa r ma n-Tre f f t z profile. " ~ 
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On transformation of thev a- plane in the £ plane 
the upper or the lower rim of "the - crescent K 1 t or any- 
one curve in its interior Becomes a slit in the £ plane. 
The straight line A^Tk^ is selected (fig. 4), that is, 
the median line K 3 ' of the crescent as slit, on the cir- 
cumference of which the tangential and normal components 
of the speed must he computed as on the Joukowski profile. 
To this end the speeds in the a plane on the reflection 
of this slit must he known, which follows from (28) at 



a = a 



s! 

4| 



sin \J/ 



1 + S. - SL 

2 



£■ jsinU'; j e 1 ^ (SI) 
8 I V J 



they are two circttlar arcs K 3 symmetrical to the u> 
axis, meeting in points A x and A 2 and having a radius 
vector of | a | 1 a, where \[/ is the circular angle in 
the a plane. After transformation the area hounded "by 



the circular arcs 
£ plane. 



K„ 



forms the second position of the 



For the speeds at the slit in the 
(3) and (9) give the flow along the u> 



d$ 



0) 



d$(a) da V 

.... ...y __ —— j i _ 



H' 



da 



• 2 



(a - a 0 )' 



I plane 
axis at : 



1 - 



t: 

a '2 



eq.ua t i ons 



(32) 



and the flow along the j 
d<P T 



di> T d$(a) T da V T 

= - i — j 1 + 

k S L 



da 



axis at : 



-1 - °~ 



- 2H 



a — a 



(33) 



o _ 



t 2 



if tan 3 = 0, a = -d , H is given hy (2)and a hy (31), 
According to (10), the source distrihution for the flow 
along the | axis follows from equation (33) at: 



' \ N> 8 , r sin ^ r d r , ■ \ , o i - 

q^-- Ty = - — V x _ ^-CcosMi. (1 - cos \|/ ) + sin 2 M/J 

a 12 

a » d 2 

+ i 5 is in \|/| cos\{/{+ 2 -- [cos ^ (l - cos w J 

+ sin 3 '4/- 3 s in 2 \J/ cos \|/] + i K Slsin^jcosi' + . . . \ (34) 

A. 1 i 1 
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and, according to equation (11), the vortex d is tr ibut i on 
of the 6ame flow direction 



(34a) 



For the flow along the n axis, equation (33) gives the 
source distribution 



(35) 



and the vortex distribution 

fn - 

s 



V - | = 1 7_ — (1 - cos \|/ ) + 2 - [1 - cosi|/- sin 2 ^ ] 



+ 3 f±_ [l _ C os\J> - 2 sin s ^+ 2 s in 2 M> c ds \J0 
as 



+ 3 - - Is In -tyl < 1 - cos * - £ s in 3 \|/ 
a 2 I I L 2 



] 



6 

+ — sin 2 \|/ (1-4 cos ^ ) + 
32 



•} 



(36) 



The transformation leaves the air flow at infinity un- 
changed. The secured singularity distributions also 
contain the circular angle \^ which is not as simply 
related to the abscissa | of £ plane as the angle 
cp in the Joukowski transformation. In this instance 
the relation is obtained by entering the coordinates of 
circle K 3 (equation (31)) in the conformal function 
(28 ) ; then 



* 2 (tan -1 ) 



a 



1 + -- 



a/k 



4 16a' 



(37) 



irt 
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cp « 2 (tan 1 ) / - 6 - a cos <p (3 8) 



For practical reasons the coordinate - — in the plane of 

a 1 

the profile is now considered as a fixed given quantity, 

f 

and \|/ , related to -- and 6 through equation (37), is 

a 1 

represented in the form 

ii « 9 + C (39) 

The circular angle cp for Joukowski profiles is given 

by equation (38); angle e is largely dependent on angle 

6, specifically € = 0 for 8=0. With (39), equations 
(37) and (38) give 



K/4 



a 



i 



- V 



«(V>"» (tan-)^i-^--i L -i-i-^2-- J - (40) 

- 1 + £ /a ' JL 1 + | /a » J 

and e(-£/a') = - e(+5/a'). The correction angle c ( -^r ■ 

V a y 

is small for the usual 6 and in close approximation may 
"be considered linear in 6 and K, respectively, as 

e (I /a') 

exemplified in figure 5, where is plotted 

against — — independently of k . 
a 1 



Expression (36) then affords with (38) and (39) the 
analytical expressions for the superposition density 
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/ * \ sin aJ/ r 

J = --.j cos \J/ (1 - cos \|/ ) + sin 2 \|/ J. 



1 - 



a 



(*0l 



+ £ (± 



1 - 2-^ — 
t a • 2 



i + e/i 



6) ■ - 



s in \|/ 



3 



1 > 

a 1 2 



Va'/I Va» /i+ 6 / a i 



sin 3 \|/ 



4V _,y 



— cos = f 4 i — ii - e 



>(4 



differing only by an additive term from the original ex- 
pression (21). They also have zero or infinity places at 
the branch point A 3 ' just as the expression, f v {% /a 1 ) , 
except that the limiting values at these places |/a' = ± 1' 



l r 



K 1 
3 J 



are now reached with the power* (1 — £ / a 1 ) p 
However, even this change in the singularity functions in- 
cluded in the correction terra is negligibly small for all 
usual angles (6 ^ 36°). 



These known expressions are supplemented by a. new 
one for the singularity distributions of the crescent: 



G?0 



= sin ^ 

i-.iL 

a 1 3 



sin y\i 



cos M» = f ( — 
\ a V 



. N 1 - 3- b — 



J 



a '2 



*If = —(a 1 — A), equation (40) becomes 

A ] i/ 2 r " ■ t 15/4 



e [-(a 1 - A)] a +2 



L 2a • _ 



I L 2a 1 J 



+ . 



18 1TACA Technical Memorandum 2ffo . 1023 



where f 5 ( -t-^N = cos cp = -L- (42a) 

\a ' / a « 

Function F 5 ^ — — J gives the source distribution of the 

cres cent 

For the singularity distributions of the Karman— 

4 

Trefftz profiles, the factor = 1 + K + , . . is sup- 

k e 

plementary to those of the Joukowski profile; it affords, 
"besides the vortex distribution caused by the thickness of 
the crescent, 

merely distributions which are small from the second de- 
gree at least. Summed up, there is obtained, besides the 
Joukowski profile singularity distributions, additional 
&is tr ibut i ons for the symmetrical Karman— Trefftz profiles: 
namely, the source distribution q 5 (| /a 1 ) and the vortex 
distribution KY 2 (|/a')» which, taken by themselves, rep- 
resent a crescent. For curved Karman-Tr ef f t z profiles the 
afore-mentioned terms are supplemented by others which are 
dependent upon edge, angle and camber. But an exact calcu- 
lation discloses all the'se terms to be small, at least from 
the second degree. 



c) Profiles with Reversed Curvature (S profiles) 

The mean line of both the Joukowski and the Karman— 
Trefftz profiles is a circular arc. But the form of the 
mean line of a general profile usually diXfers from a cir- 
cular arc for reasons of lesser center of pressure travel 
by angle-of-at tack changes. This is achieved by an S— 
shaped mean line. The deformation of the c ir _c_ular- = ar cl ike 
mean line requires an additional vortex distribution on 
substitution of the profiles by singularities. 

According to Betz and Keune (reference 3) a circle 
Kg (fig. 4) is transformed into a profile with finite 
angle at the trailing edge and with reversed curvature by 
means /of conformal function 
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+ a 



a — a 



i-cr + a 



a 



(45) 



To secure the additive vortex dis tr ibut i on a circle K x 
is merely transformed into an infinitely thin S— shaped 
mean line (fig. 6) (that is, put 6 = 0, k = 2), "because 
t.h_e different singularity distributions are known to be 
additive_in f ir s t_ appr oximat i on . Furthermore, the cri- 
terion jj> f<>r the size of" the reversed curvature is 
presumed to be small, so that terms which are quadratic 
in n are negligibly small. 

Development of the conformal function for small n 
and k = 2 affords 



i = i 'a ♦ il ♦ 1 a i °- fl - il V 



It splits in the known Joukowski transformation and an 
additional conformal function. 

The reflection of circle K 1 around 0 (fig. 6), 
a = a e ic P. is the slit, the mean line the £ plane with 
the coordinates 

£ M- . 3 T) (J. 

— = cos cp — — sin cp; — = — — sin 3 cp cos cp 

a 1 2 a ' 2 



Hence 



a ' 2 a • 



1 - — i (47) 



is the equation of the S-shaped mean line and 

tan v = L2L = - t (\ _ 3 i--^ (48) 
d | 2 V a 12 y 

the direction of the tangents of the mean line to the % 
axis. Then the tangential and normal components of the 
velocities at the slit follow at 
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v t = u - | (1-3 cos 2 cp) v 



v n = ^ (1-3 cos 2 cp) 



u + v \ 



(49) 



and equations (8) and (9) finally afford with (10) and 
(ll) the vortex distributions for the S— shaped mean line 



a' J 



x 4 y 



3 m- s in cp cos cp 



and 



V^V L y 4 X J sin cp 



> 



(50) 



In these equations (50) the fact that the conformal func- 
tions (45) and (46), respectively, changes the air— flow 
direction at infinity has been borne in mind, giving 



Vos cos a' = = Vjo cos 



Vqo sin 



a' = V T = 7a, sin (^a. - ± ) Z V 



4> Vx+ 4 Y y 



4 x 



(51) 



this change being interpreted as a rotation of the co- 
ordinate axes in the t plane. 

Then the introduction of the analytical expressions 
(21) and the limitation to terms of the first rank af- 
fords, besides the known distribution for the plate. 



( e 



2 ; = 2 V f s . — , 



the vortex distribution for the reversed curvature 



/ t \ r- - - — 

Y 4 J = - 3 H 7 X sin cp cos cp - - V x 



1 — cos cp 
sin cp 



= - 3 H V T |f 



a'/' 



+ I f 



(52) 
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the first part of which was given "by Birnbaum (reference 
=>!=)-; while the second part is" dae ,te> the rotation of the 
axes in flow direction. 

The singularity distributions thus secured f or a 
profile with finite edge angle and reversed curvature in 
first approximation, then, merely need to he combined in 
a certain manner to afford the field of flow of this pro- 
file . 

If the v_ortex distribution of .a sele cted mean line 
is to be established, it can be visualized approximately 
by the sum of circular arc (equation (6 )) and reverse 
curvature (equation (47) J. The coordinates of this ap- 
proximation can be written in the form 



JL = ( tan P - £ JL 



\ 



2a'/ 



1 - 



(53) 



Then the directions of the tangents with the | axis 
are given by 



dT) 



tan v = — = — 2 

d£ a 



L tan P - S [1 - 3 i~ ^ 
2 V a' 3 J 



(53a) 



and 
tan 
A ' 



with it the tangent in point 
u 1 = — 2tan{3 + u- 
=( I = _ a ' ) by tan V 



and the tang 
2=3 tan P + 



A x • (£ = + a ' ) by 
ent in point 



The procedure 
tan P and \x can 
(reference 1), who 



in solving the profile constants 
therefore be the same as Birnbaum's 
derived them from the tangents at 



points 



A i 1 



and A 3 1 . But, since the tangents are hard 



to establish accurately at these points, it is more prac- 
tical to define the curve of the third degree of the mean 
line from the ordinates. With ri x denoting the ordinate 

i 



at 



+ — — 

2 



and ri 2 the ordinate at \ = 



a 
IT 



(fig. 



14), equation (53) gives 



tan P = 



p, = 



2 / Tl2 + 



2Li 



£/' Ha _ ILl\ 
3 a ' a '7 



(54) 
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Shis defines the profile con- 
stants for the ease where the mea 
line is a curve of the third de- 
gree. Whether it is represented 
accurately enough "by this curve 
can be checked by the degree of 
accuracy vith which the ordinate 
value for t ■ 0 satisfies the 
condition fi 0 /a* = tan ^. 

2. Flow in Plane of Section 

a) Substitution of Profile 

by Singularities 

The flow past a profile in 
parallel flow can be achieved by 
a certain continuous distribution 
of sources and vortices on the 
mea-n line of this profile. The 
flow produced by these singulari- 
ties must be superposed on the 
parallel flow so that the closed 
streamline formed by these two 
flows agrees with the profile 
contour. The form of these sin- 
gularity distributions was se- 
cured in the preceding chapter 
for profiles of certain families. 
If these distributions are to be 
used for the substitution of 
other profiles, the extent of the 
coincidence of the closed stream- 
line with the contour must be 
checked. And this introduces the 
stream functions of the distribu- 
tions . 

The complex flow potential 
of a source element q ( £) ds in 
a point x is given by 

and for a vortex element Y(f;) ds 
by 

*.(2)=-^- In (z-i)ds 

where i is the locus of the 
element (fig. 7). Integration 



with respect to all singularity 
.elements along the arc length S 
of the mean line affords for the 
profile in parallel flow the 
stream function 

<P (z) = V m z e -,a +2^J' ? (f) In (z -Qds + 
s 

+ -^jvtt)ln(z-5)ds, (55) 
s 

with Vqo denoting the flow ve- 
locity and a the angle of 
attack respecting the positive 
x axis. The integral equation 
for the stream function 

S>(z) =y VooCOsa— x V^sm <x+ ~Jq ($) arc tg ds 

s 

+ JV (f) In i (*-«)» +• (V -W 2 ds, (55a) 

s 

is now simplified so that the 
integrals can be numerically 
evaluated independently of the 
form of the mean line. Next, 
we introduce (fig. 8) 

y » y 8 ♦ y 

whereby y s , the ordinate of 
the mean line, can, according 
to (53) s be written in form 

S~(*'-f ?)(•-£) 

Likewise 

M*'-f *)(•-*) 

and the arc element ds itself 
is given according to (53a) by 

ds^dtf^f^ 

the small quadratic terms being 
discounted . 

Herewith the integrals be- 
come +0 , 

JV (f ) In V (*-£)«+(»— nY ds^Jy (£) In V ( x - S )*+y* d c 
s ^ 

a , igP J Mil {x _ s) t + yt 
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+»' 4 (*-*)(*+*) 



aad^ 



— a' 



(*-f)» + ? 



(66 a) 



arc tg 



arc tg 



(66b) 



— a' 

In addition to the underscored 
integrals independent of the 
form of the mean line, we obtain 
cumulative integrals multiplied 
by consistently small values 
tan 0 and |x and, together with 
the singularity distributions, 
they are therefore small, at 
least of the second degree. Hence, 
they have no further significance. 
We the-rel'ore dispose, similarly 
to Birnbaum (reference 1) ana 
Glauert (reference 2), the singu- 
larities on the chord of the mean 
line, and so obtain 

+*' 

!P (z) = y V m cos a — x Fco sin a + ^ Jq (f ) arc tg x ^_^ d£+ 



section must, of course, 
carry the same identification. 
But, since this stagnation 
point is not generally known, 
some other method must be 
applied* The flow is there- 
fore divided again into its 
components along the axes of 
the coordinates and the con- 
stants are so defined that, 
for the flow in the x direc— 
ticn, the rear stagnation 
A}' and for the flow in the 
y direct i on, the profile point 
on the j axU (x - 0, y«y"i) 
lies on the streamline » 0. 
Then stream function (57) 



gires 

fx y 



a'V x 



+ 6 [ W *& ^)-*4<l;0)]} + fa (58) 

for the flow along the x axis* 



and 



(57) 



for the stream function. 



According to the fore- 
going, the singularities ex- 
tend from A x * to. A B '. 
Putting = const affords 
an integral equation indicat- 
ing^ the field of "all stream— 
lines in the plane of the 
section, one of which is the 
clos.e.d_s treamline , that must 
agree with the section con- 
tour. Since all streamlines 
are defined merely up to one 
constant it is decided that 
this streamline shall obtain 
the ralue - 0. Then the 
streamline of the parallel 
flow which goes in the for- 
ward stagnation point of the 



+ 2 4Nv= • - (59) 

for the flow along the y axis, 
with abbreviation 



— a' 

for the source— sink integrals** 
and 

♦The — denotes the correction, 
according to formula (42), 

♦♦formula (60) is evolved from 
the relation 

+ (tan~ 1 )-i r «TT-(tan~ 1 )^' 
x-S y 

whereby, because of (23) 



■n a 



-a » 
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n a 
— a' 



for thd vortex integrals, and 
g x and gy denoting the quad- 
ratic small terms. ?or the 
profile itself, we get ilf _ = 0 

*y - 0 















da; 


bx 



(63) 



(62) 



If the four constants — , 6, 

a 

tan 3, and |x are correctly- 
defined the contour of gener- 
alized Karman-Tref f tz profile 
(reference 3) should agree ap- 
proximately with the stream- 
line from equations (58) and 
(59) , if necessary with allow- 
ance for the quadratic terms. 
The agreement can he checked 
numerically hy means of the 
appended graphs which give the 
values of the flow function 
for the source distr ihutions 
<iv (equation (60)) and for 
vortex distributions W v 
(equation (61)) plotted against 

y 

the ordinate — for various 
a * 

abscissas ~. 



b) Field of Velocity in the 

Vicinity of a Profile 

As a consequence of the 
resolution of the singularity 
distributions which transform 
a profile to a streamline in a 
parallel flow, the velocity 
field around this section is 
established without difficulty. 
The components of the speed are 
given by 



— it, 



b x 



and, according to (62), 



In this case also the effect of 
each singularity distribution 
is given separately, in order 
to preserve the lucid aspect of 
the effect of the individual 
profile constants. 

Partial integration of (58) 
and (59) gives the four integrals 



bQ v _ 



r— £ 



b y/a' 



bxja' 



(x — tf + y 



--i{*MM*)MM*)]} 



+ 0' 



l nS h W) y* d * 



+ a' 



by/a' 



and 

b W v 



l $ fr [i')( X - 



71 

— a' 



?) 2 + y 



a t/a' ~ I j f '[l)lx — W + p de 

which divide into the two inte- 
grals 

— a' 

and 



The singularity distributions 



J 



can be written in the 



form 



Therefore the resolution of the 
one general integral 



(64) 



(65, 
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— a' 



' + f (x— £)' + 

-fco' 



/i\»-»r df ehange. sign, the root being 

= — — ill— . tff =-, . . (66) posit ire for x < 0 and nega- 



is sufficient, if n is a whole 
positive digit (n » 0, 1, 2, 3). 
After the substitution t, - x » 1 



a'— X 
-<a'^x) 



d to 



y a '»— (» + *)» 



tire for x > 0. 

Denoting the integral J n 
for the individual values of n 
with J 0 , J lt Jz, and J 3 
while keeping fyU/a 1 ) in mind 
equation (65) gives 



Development of the integrand of 
this integral for n = 3 in 
powers of w gives the integral 
in the form 

(a')"" 1 J = b 1 L 1 + b t L t + b s L 3 + L t . . . 

where 



= f J» „ + , 

J ya'»— (* + «)» 
(a'-Fx) 

a'— x 
( wdw 

J 1*'* — {* + *)■ 
-(a' + x 1 

= ° r = + + ^ , 

J ya'i-(* + ir)» \ 2 / 

-(a'+x) 



■Pi 
-Pi 



(67) 



and 



Ml) 



Mi): 

'[''(?)] '■ + '•! 



= J x + 2 J t 

= A 



p '[Hir)\ = J ' + J * J 



and 



(68) 



which, entered in (64), finally 
yields the dimens ionless speed 
factors of the vortex distribu- 
tions, as follows: 



Along the x axis 



while the solution of integral 
leads to 

J w s + V* f«' 2 — (*+«»)» a' 
-(o' + x) 

where 



aw, 


y 


ay/a' 


a' 




y 


dy/a' 


~~ a' 




_ y 



and 



yi[»-(>-*+s)i 



aw , 

a a/i 



are dimens ionless values. It 
should he noted that 



1 + 2 * fll + ( 1 _£+JL) fli ] 

-£(7*-»^ + £(»7*-*)] 
and along the y axis 

&='-('-5h + S°- 

&-7-('-^)°-+S(^ D -- D .) 

The dimens ionless speed values 
of the source distribution are 
obtained from those of the vor- 
tex distribution with (22) and 
(64) . The value for the com- 
ponent along the x axis is 
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*y /a' dx/a' 

and that along the 
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(74) 



-2- !, (equation (42)), the 



d 2/a' " d y/a' 8 y/a' 

How there is lacking only 
the speeds of the source— sink 
distribution of the crescent 

"\ 

a- J' 

flow function of which occupies 
a special position among the 
integrals. Its exact evalua- 
tion requires a numerical proc- 
ess, since f e (|/a*) is given 
as such (fig. 9)*. But, dis- 
counting the region in proxim- 
ity of the branch points A x ' 
and A 2 ' , the analytically 
given form f 5 (|/a») (equa- 
tion (42a)), which is wrong only 
in the vicinity previously men- 
tioned, can be used (shown 
dashes in fig. 9). Then 

by/«' n a' J (x — f)»+y 2 * * 



by 



and 



+ a' 



±±r 

it a' J tx — 



d x/a' ' 7i a' J (x 
— a' 

Solution of 



(75) 



f) z + S ! 

the two inte- 



grals 



J, 



{x-W + y 1 

-a' 

and 



. . (x- <0« + »» , *. ^ te 



1, 



a' J + i/ 2 ~ + «' (* + a') 2 + V 



+ 



■ —_?/'_ 
a' i/ 



arc tg - 



2 a' y 



gives for the source distribu- 
tion of the crescent the func— 
t i ons 



♦The values fT 



( s 



can 



in 



; 

good approximation, be considered 
linear in 6 » K n and are 
therefore plotted for 6 ■ 1.0. 



0 y/a' 



<) x) 



a' ji I 2 a' 



^ 2 a' 1 " + p 

— y ; arc tg - 



In 



(x — a'Y + Xf 
(x+aV + V- 



2 a ' y___ 
— a' 2 + y 2 



2 a' )/ 



+ — arc tg — r .— 

a — a 2 +y 



> (70) 



Herewith the dimensionless 
speed factors of all the singu- 
larity distributions are secured 
which, in first approximation", 
give the velocity field in the 
vicinity of a profile. The sec- 
ond part gives the velocity fields 
of every single singularity dist- 
ribution in diagrams which euable 
easy reading of the values for 
each point o The individual values 
obtained by differentiation of 
equations (58) and (59) ai.a al_,„- 
ance for (63), read as follows: 



== uos a 



a \a a 

8Q.(4.4) »w 3 ( J , T ',) 

»»■.(;■•*) 



0 y/a' 

I) 



+ sin a 



3 y/a' 



C — 

i) y/a' 

aw/*; 7/ ) 
\a' a') 

9 y/a' 

a w,(4;4 

3 \ ra" ra' 



I a&s_ 
^ a y/a' 



a £> y/a' 

speed_ in y axis direction 



"a y/a' 



4* -4) 

la a'/ 



= sin a 



I -(1 + 2- + xi il-J II 



- cos a 




+ 2tg/» i!L_?J 
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Here also the terms with g v and g are quadratic 
small (equations (82 ) and (83)) in the appendix. 

Because of the cited omissions the formulas secured 
here are not applicable at the profile edge or its imme- 
d ia t_e vie ini ty , wh i c h i~ iiior e over ', w ouTcf exceed t he pur- 
poses of the present work. 

3. Appendix. The Quadratic Members of the Development 

Since the singularity distributions were obtained 
by means of a development according to profile constants, 
the linear terms of the development will ordinarily be 
insufficient in closest proximity of the profile contour. 
This fact is of importance on_ly when attempting to check 
the accuracy of a predetermined profile by the given 
distributions. But for predicting the velocity field in 
the neighborhood of a profile the linear terms should 
prove sufficient, whence in this case the quadratic terms 
give in error estimation. 

Quoting all the quadratic terms secured would com- 
plicate the representation" and entail hardly worth while 
paper work in the error estimation and the check on the 
profile contour. Hence the various nonessential terms 
are disregarded. These are the correction terms in equa- 
tion (41), which secure new analytical expressions for the 
superposition density. They are themselves small as quad- 
ratic terms and produce minor, negligible changes in the 
distributions of the Joukowski profile due to the edge 
angle. The correction term in equation (42) alone remains, 
in order that this distribution portray the true conditions 
of the crescent of the Kar man— Tr ef f t z profile. The addi- 
tional integrals in (56a) and (56b) are discounted, as the 
new integral forms have only a slight effect on the results. 
All these terms can be disregarded for the very reason 
that the general profiles, for which ultimately the singu- 
larity distributions are used, mostly deviate from Karman— 
Trefftz profiles from which the distributions have been 
obtained. 

There remain then the quadratic terms of the singu- 
larity distributions of (16) to (19) and (34) to (36). 
Of these only the distributions multiplied by V ^ = 
Vco cos a are employed; the quadratic values multiplied 
"by Vy = Yoa sin a are likewise ignored, as the angle 
of attack a is consistently presumed to be small. 
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Furthermore, while the wing chord is largely responsible 
for the flow in y— axis direction, that is, the vortex 
distribution of the flat plate (equation (24)) the thick- 
ness and camber of the profile are of secondary importance 
In consequence, only the singularity distributions are 
written that appear for the flow in x axis direction, as 
they alone contribute a perceptible share to the results. 
They are the source distributions 



/ 



d 

a 



2 V x - ( 2 - + K , f , ( 

x a V a J Va ' 



\ 



' Jt ^ 
J 



\ 



= + 12 7 f 
x a 3 



a 1 J 



> 



(79) 



- (*){ t 

4 6 ' " 



i ]= <- — V 6 k f . C-^-^ 

\ a •/ 2 x 5 Va '/ 



of which the source distribution f. 



, ; — r did not appear 

4 V a 7 

previously but will gain significance for th<. design of 
new profile forms because of its property of shifting the 
rirea of maximum thickness of a profile more toward profile 

(a)/ % \ 

; — ■ and 

Q. K ( — t ,' afford a correction to the contour of the 

\ n ■ J 

Joukowski profile and the crescent, r e s pe c t i ve ly , Because 
of the thickness effect on the cumber the- vortex distribu- 
tions 



center (II, 5). Distributions 



(s) / l \ 



d ( I N 

- tan 0 f 1 & 



Va .' y 



(2)/ j_\ 

Va-y 



Y_ , " i = + 4 T x ^2 | + K J tan & f 3 ^ i- j v, 



12 V x a tan 0 f. 



a 



(80) 



are additive to the corrections. Without these quadratic 
small terms (80) the profile on the side facing the curva- 
ture radius of the mean line would be thicker than on the 
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other (fig. 21): that is, the profile would apparently 
.hay e a s ma ller camber. 

In this manner the correction term for the stream- 
line of the profile by V x flow is obtained: 



• 8 r , "'[ , '(r i ff)"" ttlilll j' (B1) 

and inserted in equation (58), while g (equation (59)) 
is "being disregarded. 

Differentiation of this equation affords the quad- 
ratic small uortions to the speeds in equations (77) and 
(78): 

5 v/a ' aV a y 5 y / a i a* S y/a ' 



&y/a« aV a J h y/ 

° H5 \ _ 1 » . . / j 1 " ■ a 

4 by/a' a ay/ 

+ 2 ( 2 i +K tanP -StllilZ -S^taaP (82) 

y by/a' a a y/t 



+ — K 6 - i tan p — 



and 
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the individual terms of which are known from the linear 
distributions (equations (?2 ) to (74)). Newly added are 

the speeds of the source distribution f .[ — i which, 

Va' y 

however, can be derived from those of the known vortex 
distribution for the S— shaped mean line, according to 
equation (64), Because 



S Q 4 f^;i-) & w/^iV, 
__\a_^aJ_y = V a ' a ' J 

d y/ a 1 d x/a 1 

and / -r T? ^ / x v N 

* ( ~;-t ; a w 4 ( - x T ,.- y 7 ; 

\a ' a 1 J \a ' a ■ / 

S x/a ' by/a 1 



(84) 



(85) 



d Sv ^ g v 
The terms 11 — and — - — are disregarded conformally 

d y/a 1 d x/a ' 

to the foregoing arguments. 



II. APPLICATION 
1. Generalities 



The solution of the singularity distributions for a 
generalized Ka'rma'n-Treff tz profile (reference 3) is de- 
scribed in part I of this report. The outstanding fact 
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is that the vortex and source distributions necessary for 
the~ sub's ti'tut ion of a profile can be combined from only 

five different distributions f« {' ~ indicating the 

V a ' J 

abscissa of one singularity. To make the calculation as 
much as possible independent of the angle of a.ttack a 

the parallel flow Voce"* 10, is divided i,nt o its component 
cos a along the x axis and its component V y = 

Voo sin a along the y axis. The first result is the vor- 
tex distributions as previously given by Birnbaum (fig. 

y 



10): the vortex distributions 2 V v f B f \\ of the flat 

V a J 



f t N 

plate, the vortex distributions 4 V_ tan P f 3 i ~ , ■ of 

N " ^ ' ' t \ 

the curved plate, the vortex distribution —3 V_M,f 4 ; — — , 

of S-profile.* These distributions are supplemented by 
two new source distributions (fig. 11) which appear only 
for the fldw~~"V x .' The one 'source distribubion 

d ~ ' f s ; 

— 2 V„ — f . ( -i— leads to a closed streamline similar 
X a V a * J i 

to the symmetrical Joukowski profile, where - is the 

parameter for the profile thickness. For the Karman— 
Trefftz profiles with finite angle 6 at the trailing 



edge the source distribution \—V_ S f 5 ; -=7 , is additive; 



X 5 \3J J 



which, considered by_.it self, is the source d i s t r ib ut i^ n_ 
of. a, symmetrical crescent" w_ith~edge. angle 6. (fig. 11).** 
Besides these singularity distributions generalized 
Karm&n-Tref f tz profiles have still other distributions 
composed of these, five basic funct ions, which must be com- 
bined with the af ore— ment ioned superpositions and with , 
the parallel flow according to a certain law if the Karman- 
Trefftz profile is to be streamline in parallel flow. 



•Constants tan P, ij. , — ., and 6 are discussed farther on. 

a 

/ t >. 

**e; — . is explained in I ,. 1 b ; for the rest, the graph 



affords a satisfactory representation of the aspect of 

/ 
/ 

V a' J 
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Then, in order to establish the effect of any one 
of the singularity distributions on the flow formation 
or the 3-peeds' at a point of the plane of the section, the 
effect of all its singularity elements at this point must 
"be summated, that is, integrated along the line on which . 
they are located. The mean line of the profile is chosen 
as the depository of all singularities. During the devel- 
opment of the integral it was found that the ordinate 

a » 

of the mean line (fig. 8) could be disregarded, and the 
singularities themselves assumed on the x axis, with a 

y e 

much simpler evaluation as a result. With — as the 

a 1 

ordinate of the mean line, the point of the profile plane 
is given by 




(66) 



with - — denoting the distance of this point from the x 
a 1 

axis (fig. 12). For starting points at greater distance 

y y 

from the profile, put — e -A— , that is, let the mean 

a » a ' 

line and the coordinate axis coincide; but the two values 
must be carefully differentiated in profile proximity. 

When in the following a method for computing the 
plane potential flow about general thick wing profiles is 
evolved, it is practical to utilize fully the data secured 
for the generalized Karman-Tr ef f t z profile, since the flow 
condition for each of these profiles is known. As the 
calculation is to deal with the flow about a general pro-' 
file- in the neighborhood of the profile but not on its 
surface, the method can be restricted to an approximation, 
that is, the given general. pr of i le is closely substituted 
by a generalized Kar man-Tr ef f t z profile, which merely 
stipulates knowing the profile constants for this similar 
Kar man— Tref f t z profile, which then is a good approximation 
for nearly all practical profile forms. The substituting 
singularities of this profile are herewith known, and it 
will not be necessary as a rule to check on the given 
profile the extent to which it is represented by the 
s ingular it ies . 

The reason for checking the profile contour (II, 3) 
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is simply in order to show how closely the contour of 
^g.ejigr.al profiles agrees with the closed streamline of the 
approximately resolved flow.' 

But the singularity distributions afford, in addition, 
still other new profile forms with certain, predetermined 
profile characteristics. They may even vary substantially 
from those of the Kar man— Tref f t z profiles if provision is 
made that the ' s ingular ity distributions characteristic 
for certain profile forms are combined in a manner other 
than that prescribed for the Kar man— Tr ef f t z profile. The 
profile contour is computed from the singularity distribu- 
tions analogous to the method employed by Fuhrmann (ref- 
erence 6) to airship hulls. By the determination of the 
velocity field of these profiles the pressure distribu- 
tion can be indicated also in good approximation.* 



The generalized Karman— Tr ef f t z profile is completely 
defined by four profile constants. But before the calcu- 
lation can be made the mean line and the axis of the 
Kar man— Tref ft z profile must be plotted first in the given 
general profile, which agrees with the mean line of the 
general profile only for symmetrical profiles. The mean 
line of the section is closely approximated by the line 
connecting the centers of all the circles which touch the 
profile contour twice from the inside (fig. 13), the cur- 
vature circle (radius p) on the nose included. The axis 
of the Kar man— Tr ef f t z profile which at the same time is 
the x axis of the plane of the section is given by two 
points A 2 » and A 3 » (fig. 13). Point A a ' , located 
at the section tip, forms the vertex of the pointed trail- 
ing edge; for sections rounded off aft** of the rear con- 
tact point of the tangent to the section lower surface, 
of the chord, with the section is chosen. Point A 2 ' is 



the mean line.*** Distance A 1 , A a l is put equal to 2a'. 



*Although -it affords the closed, streamline for the flow in 
x direction exact , the vortex distributions for the flow 
in y direction can only be approximately determined. 

**Strictly speaking, it should be constructed just like 
point Aj 1 at the-nose. 

***Herewith the position of point A 3 ' is accurately 
enough defined; although the plotting of the curvature 
circle is afflicted with a certain error. 



2 . 



Determination of the Profile Constants 



located at half 




from the nose on 
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The geometric angle of attack a g does not, in general, 

agree with the angle of attack a which must "be measured 
from the x axis.. With a B denoting the angle between x • 
axis and profile chord (a g = 0) (fig. 13): 

a = a g + a s (87) 

The ordinate axis, axis y, of the profile is the median 
▼ertical on the length A x ' A3 1 . 

Quantities tan 0 and |x giving the camber and the 
reversed curvature are obtained from the ordinates of the 
mean line and not from the tangents (as Birnbaum did) 
which cannot he accurately plotted. For generalized 
Ka'rman— Tref f 1 2 profiles the mean line is, in good approxi- 
mation, a curve of the third degree (equation (53)) through 
points A 1 'Ag', with which the general profile must coin- 
cide in two additional points. Choosing the points at 

-X- a -0.5 with the ordinate ~£ and at ~ = 0.5 with 
a ' y a ' a ' 

the ordinate -4 (fig. 14) affords* 
a ' 

1. The index for the camber 

.„,.£'!•♦ If\ (68) 
3 ' a ' / 

2, The index for the reversed curvature 

, = */■ Zl - l±) (89) 
3 V a • a 1 J 

This secures the quantities which define the intensity 
of the vortex distributions for the cur ved_ S— shaped mean 
1 ine . 

The two remaining profile constants of the general- 
ized Karman— Tref f tz profile must be computed from the 
thickness curve of the profile; whereas a Joukowski pro- 
file requires only a section constant d/a for the thick- 
ness, the finite trailing edge angle 8 on the Karman- 
Trefftz profiles affects the thickness at every point, 
particularly, a profile of finite thickness, the crescent 



*See also: equation (54) and the context. 
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d 

is left over for — = 0. Therefore, the thickness of a 

a 

Karman-Tref f tz profile is defined by two criterions, 5 

and — , . the magntiude of which is not ascertainable from 
a 

the plotted profile without specified assumptions. If 5 

d / 

were given the value of — , (denoted with by Betz — 

a a 

Keune) could- be deduced from the thickness at profile cen- 
ter by the Betz— Keune method (reference 3, equation (2)); 
but it is found that 6 cannot be accurately enough de- 
fined from the angle at the trailing edge, because the ef- 
fect of the constant — on the profile thickness starts 

a 

too close to the trailing edge, as is readily seen from the 

Karman-Tref ftz profiles.* Therefore, 8 and - had to 

be secured in a different way. Since the thickness remains 
approximately independent of the camber and the S curvature 
(see equation (2) of reference (3)) and since the symmetri- 
cal Kar man— Tr ef f t z profile acts, in good approximation, as 
if the ordinates of the crescent with the edge angle 6 
and the ordinates of the Joukowski profile with the cri- 

' d 

terion — were additive,** the profile constants can be 
expressed by approximate formulas. If profile criterions 

are employed again at — • = +0.5, it then is expedient to 

a ' n 

D , x 

denote , as thickness at — = 0,5 and — as thick- 
s' ' a ' a ' 

ness at — « -0.5 (fig. 14). 
a ! 

Then 8 follows at 

- * -( ) |73 + ( I (90) 

a 3 ^a ' a » ' L. V a ' a 1 /J 

and 6, in radians, at 



*Angle 8 at the trailing edge is read off as much as 80 
percent too great. 

+ *For the symmetric Joukowski profile and the crescent the 
ordinates in each point are easily secured, but that for 
the Karman-Tref ft z profile is very complicated. 
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6 ~-i 3 (91) 

3 V a 1 a ' / 

Although these formulas are not exact, equation (90) 
proved very satisfactory in a check on a large number of 
Karman-Tref f tz profiles, and equation (91) agrees with the 
correct values up to K 1°. Since equation (90) is inde- 
pendent of 8, an error in 6 is of little consequence 
on the subsequent results. With these constants the aero- 
dynamic characteristics of generalized Karman— Tref f t z pr o— 
f'TTes themselves are secured. These relations also hold, 
appr oxima£.e_iy , _f?r general profiles with the same con—" 
stants. The theoretical lift" is 

c„ s - 2 tt / 1 + - - + sin fa + 0 - i (92 ) 

a t/a' \ 2 Tt a y V. 4 * j < 

where ---- (\ + — — + — is the theoretical lift gra- 
t \ 2 tt a J 

dient that is always greater than the true value, because 
the friction is not taken into account, and a 0 = 

- • i \x is the zero lift'angle of this profile, 

measured from the x axis. The moment coeff icient, referred 
to the centroid F of the profile, is (see reference 3, 
p. 341) 



i 1 6^ 

TT 1 — — — 



AJL— 3-Jli- S in 2 - | (93) 

8 



t 2 /a' 2 



The profile has a fixed center of pressure, when Hj^— £.* 

3 

3. Substitution of a Given Profile by Singularities 

In a check on the extent to which the generalized 
Karman-Tref ftz profile obtained by the computed constants 
agrees with the given general profile, the study can be 
limited to the profile in parallel flow along the x axis. 
The parallel flow' along the y axis is essentially affected 
by the vortex distribution of -the fla't plate into which 
the profile chord enters. Against this effect the dis- 



*See reference 3, equation (48), where, however, the 
squared terms- are taken -into account. 
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tr i"but ions due to profile thicknees or to the mean line 
are inferior, especially since V y always remains smaller 
than V^. (Vy V x ) ; hen c e the approximation of the profile 

"by one of the Karman-Tref f tz type should always remain 
accurately enough in this direction of flow. 

For the flow along the x axis the streamline equa- 
tion (58) — if the closed streamline is given the value 

— S_ a o * - affords 
a'V x 

. a , a "\a''a-/ 2 S\,a' a' 7 

+ 2 tan ; - »,(U0)j- S ^L;Z. V ».Cl|0 ) 

+ 6|W 4 / JL..JL N j- W 4 (l;0) |> + g x (94) 

Here Q v ( — j and W v ( ; JL N j are, according to 
\a a y ^ a a J 

equations (60) and (6l),' the abbreviations of the inte- 
grals of the flow function of the respective source— sink 

and vortex dis tr ibut i ons f,| The subscript v 

/ t \ Va' J V 

of fy ! — - ) is always in agreement with that of Q, and 
\a J v 

Wy. ' These integrals have been evaluated for all neces— 
xv 

sary values — and — — and the results plotted in fig— 
a 1 a ' _ 

v x 
ures 15 to 19 against — — for fixed abscissas — and 

a 1 a ' 

tabulated in tables 1 to 5. Allowance for the squared 
term g x 



*Provided this streamline in the rear stagnation point 

of the profile A i'( ~ = +1.0; = 0^ reaches the 

value = 0. This is achieved with the aid of the 

a'V x 

constants W B (1;0), tf 3 (1;0), and W 4 (l;0). 
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3 - /xyN^d r / x y \ 

+ 2 ^2 * + k S ) tan P jw^JL.Z.^ _ W,(l;0) 



i ua.ii h w i — — , — — i 

ay l 3 Va''a'y 

/ x y \ 



6 - tan P |W 4 < - W 4 (l;0) 

d i_ \a 1 a '/ 



(94a) 



/ X y \ 

adds the integral Q, . — ; — ) of the source-sink dis- 

V a' a' / 

tribution f<i(|/a') the values of which are shown in fig- 
ure 20 and table 6. The squared term is usually needed 
in a check of the profile contour. 

The streamline simulated by the singularities is de- 
fined by approximation on the assumption that the given 
profile actually has become the streamline through the 

computed profile constants tan 0, n , — , and 6 (equa— 

a 

tions (88) to (91)). In this instance the right- and left- 
hand sides of equation (94) would have to agree. Then 
the successive entry of the coordinates of individual 
points of the profile contour in the right— hand side of 
(94) with due allowance for (86) — the values of the in— 

tegrals Q ( ~; — l and V I — ; — j being read from 
v \a 1 a ' / v \a 1 a ' / 

figures 15 to 20, and the calculation proceeding in this 
manner with (94) — ultimately affords a value, which shall 

/ y \ 

be expressed with — f — - } . It affords already a very 



a 



good approximation for the ordinate ( — ) of the 

\ a ' 

closed streamline -fy x - 0. 

Since all singularity distributions were directly 
derived for the contour of a generalized £ar man— Tr ef f t z 
profile, the form of the given profiles is approximately 



NACA Technical Memorandum No. 1023 



39 



always complied with "by the singularities. Hence, no 
pecular it ies are anticipated at any profile points, and 
it is sufficient to compute the ordinates of the stream- 
line in a few points. 



The method is first checked on three Karman— Tref ftz 
profiles (fig. 21) constructed according to the Betz- 
Keune method (reference 3) with the profile constants 
given in figure 21.* The constants are computed accord- 
ing to equations (88) to (91) and the values posted in 
equation (94). The result is 



Z JxflA-juLaLLLa : 
(b** = 0.138 X ) 

V t J 



— = 0.272; -- = 



a 



d o 
- = 0.089; 6=9 

a 



i.£jQj3nd._pr^i£.i2A : 
m 0 .120j 



\ t 



=-S = 0.234; -- = 0.112; 



- = 0.075; 6=8° 

a 



0.130; -- = -- = 0.0 
a' a' 

tan 3 = 0 ; = 0 



; . Z2 a ll = 0.038 



tan P = 0.051; »a = 0 



Third profile : 
(?J£* = 0^195) 



— = 0.3 98; -- 
a 1 a' 



0. 130; — = — = 0. 066 



- = 0.15; 

a 



6 = 8 



tan 0 = O.088; \t = 0 



Comparison of these values with those in figure 21 shows 
a slight error in edge angle, which, however, is scarcely 
noticeable in the calculation since the effect of thick- 
ness d/a is much more decisive for the profile contour. 
Figure 21 indicates that the second approximation (squared 
term g x accounted for) is a very good check, even for 
great thickness, of the extent to which a given profile 
agrees with the Karman— Tref ft z profile. 



As an example for general profiles, three 



files are selected (reference 7) and 
constants are computed (see fig. 22): 



first the 



1TA CA pr 0- 
-pr of i le 



d 

a 



a 



and tan 3 



00^ 



in reference 3 
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NACA 2412: 



» 



-i = 0.241; ~ » 0.125; ~ = 0.0323; ~ = 0.0246 

a 1 a 1 a 1 a 1 

- » 0.072; 5 = 10° tan 0 = 0.038; n = 0.021 

a » 



"NACA 2712: 



-2. = 0.240; = 0.129; = 0.0242; ^ = 0.384 



a « 



- = 0.068; S = 11° tan P = 0.042; u = -0.038 



C- 



a 

NACA 4418: 

-5 = 0.386; ~ = 0.190; — = 0.0667; ~ = 0.0525 
a 1 a a 1 a ' 

d o 

- = 0.109; 8 = 16 tan 0 = 0 . 0"8 ; u = 0.038. 

a 

To explain the calculation method, the calculation of 
y \ 

- j is carried out for a point of the NACA 2412 air- 
foil section. The profile point selected is = — 0.5; 
_ * a ' 

= 0. 153 . Then (86 ) 

a ' 

-2- = 0.153 - 0.032 = 0.121 and there is obtained 
a * 

^!(-0.5; 0.121) = -1.11 W s (-0.5; 0.121) = -0.308 

Q 5 (-0.5; 0.121) = -0.298 W 3 (-0 . 5 ; 0.121) = -0.033 

[Q 4 (-0.5; 0.121) = -0.195] ¥ 4 (-0.5; 0.121) = 0.167 

V a (1.0;0) = 1.0; W 3 (l.0;0) = 0.25; W 4 (l.0;0) = -0.166 

. , / y ^ 

wherewith equation (94) gives: — ! — — , = —0.131 + g . 

The accord with the ordinate of the profile point is still 
unsatisfactory. With the quadratic term (equation (94a)) 

g = -0.018 "becomes ( ~ ) = 0.149 and affords a sat — 
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isfactory agreement "between profile point and streamline 
of~ the 'related generalized Kar man— Tr ef f t z ' pr of ile . The 
same procedure is followed at the other points of the 
section contour. The results, plotted in figures 21 and 
22 manifest satisfactory agreement at all points. 



4. The Velocity Field of a Given Profile 
The singularity distributions on the mean line of a 
profile induce at each point P ( ^~ — ), (fig. 23), a 



i * 



\a ' a ' J 



speed w = w e" 1 "'", the slope of which toward the positive 
x axis * and the absolute amount Jw| = w. The amount 
and the angle follow from the components parallel to the 
coordinate axes as follows: 



w 



VV^ / VVco J To, V 2 



J 



and 



= (tan -1 ) 



•u u 



(95) 



(96) 



The speed components — and 

Too 



are obtained by 
differentiation of the equation of the flow function 



Too 



a 



Tq., with respect to 



and 



being composed of 



the differential quotients of the flow function ^ x /a ' Too 

obtained for the flow along the x axis, and the flow func- 
tion ^ y / a ' v oo along the y axis; it is: 









f 


u 


cos a 




a' T 

Stan a 


v« = 




by/a' p 


y/a' 


_y_ _ 


— c os a 


" *x 
b ?— b 

__^i_Y.2? + __ 


_ *jr__ 

±l_Zs°tan a 


Too 




b x/a' b 


x/a ' 


speeds 


u , 
y— and 


--, being 


thus known 



(97) 



00 
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tions (58) and (59), can be defined by the method given 
in the tabulation. Bach source and vortex distribution 



) and tan 0, n, — 
a 1 J a 

utes its share to the speeds (made diinens ionless with V^) 



composed of fy^-^J and tan 0, n, — and 6, contrib— 



at a point P ( — — - ), for the source distributions 
\a ' a 1 J 

!>■ 

•■SuLZ along the x axis and 



a 



i 



6 % ("^) 
Va ' a ' J 



along the y axis 



d x/a' 

Correspondingly the vortex distributions are 

V_a^_aV/_ and Va' a ' 7 

5 y/a 1 d x/a 1 

Thus it affords, according to the tabulation, for example: 

a ' Too 1 + d + 

d y/a 1 a o y/a 



• • • • 



The velocity field for each of the five singularity dis- 
tribution is shown in figures 34 to 28, the lines of con- 
stant speed in x and y direction being connected in 
a system, from which the speed portions of a selected 
singularity distribution for each point can be read off.* 
The calculation of the correction terms is usually un- 
necessary, since they merely afford an error estimation 
and are given in (82) and (83) simply for the sake of 
completeness. One exception is the quadratic term 

, which essentially affects the direction of the 

d x/a' 

*Beeause of the necessary additive term (equation (42)), 

the speeds in the vortex distribution f ( -L ) (fig. 25), 
were obtained by numerical integration. 
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o ^ x /a ' V Q 

speed. . This is due to the fact that the speed 



b x/a' 

to "be multiplied "by cos a- begins with terms already 

b g x 

small of the first order. Here is of the same 

order of magnitude as the linear terms of tan a 



b x/a ' 

As an example, take the velocity field of the 
G-ffttingen 624 airfoil section (reference 8) (see fig. 29). 
Computing the profile constants by equations (88) to (91) 
af f ords , with 

-2 = 0.319; £i = 0.160 
a » a' 

7s y i 

-7 a 0.102; -- r = 0.062 
a ' a ' 

the values 

- = 0.10; 8 a 12°; K = - = 0.067; tan P = 0.109; p. = 0.107 
a 

hence for the zero lift angle: — a Q = P — i p.; a Q = —4.65°. 

Since the profile chord slopes at angle <x s = 2.2° toward 

the x axis, it is found (equation (87)) that the zero 

lift angle a 0 = -4 . 7°-2 . 2?-6 . 9° is in complete agree- 
ment with the experimental value. 

The speeds at a point of the plane of the section 

x y 

with coordinates — = —0.5 and — = —0,3 are to be 

a ' ' a 1 

computed. Figures 24 to- 28 give 

d Q, o q 5 d W a 

= 1.1. -- = 0.23 -- = -1.52 



o y/a 1 b y/a ' b y/ a 1 



b W 3 b W 4 

1 -= 0.63 ■ ?= 0.21 

by/a 1 , by/a 1 



b Qi S (l 5 b W 



2 



= 0.25 = 0.-31 = 0.42 

b x/a ' b x/a ' b x/a ' 

b W 3 b W 4 

= 0.34 = -0.15 

b x/a ' b x/a « 
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for which the tabulation accordingly shows the following: 
~ = cos a [1.006 + -1.806 tan a] 

CO 

V ^ g x 

=- = —cos a [-0.005 — 0.397 tan a] — T ~. — cos a 

J d x/a« 



and 

J S 
b x/< 



x =-0.018 



* 



The only quadratic term to be considered in the caleula— 

tion is » — . For the mean line a^fcone , which up to now 

b x/a» 

indicated the substitution of the profile according to 
Birnbaum, the same calculation method yields 

«- = cos a [0.871 + -1.53 tan a] 

0= 

▼ h %x 
— = -e-cos a [—0.062 — 0.58 tan a] — cos a 

Vco ox/a' 
and 

d x/a» 

With a = 3.4° for the angle of attack which, according 
to (92), gives the lift c & = 1.0, and the profile thick- 
ness taken into consideration, it is found that 

~- * 0.897; ~ = 0.02 9 + 0.018 = 0.047 

. ^CO CO 

and for the mean line alone, according to Birnbaum's more 
convenient method in this instance 

~- * 0.779; ~- = 0.097 

CO CO _ 

y 

This calculation is continued for other ordinates at 

a 1 

point — = —0.5, with added choice of angle of attack 



a 
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a = -4.65°(c a = 0) and a = -0.7°(c a = 0.5), The re- 
sults are shown plotted against y/a ' in figure 30; 

whereas figures 31 and 32 shows other sections for = 

a 1 

constant. In all these graphs Birnbaum's approximation 
is indicated with dashes. The solid curves give the 
values with profile thickness allowed for, including the 

9 g 

quadratic term tn the dash— dot curves this 

ox/a' 

term has been discounted. In figures 30 and 31 the ve- 
locities were not computed up to the profile contour 
"because the approximation is no longer applicable there. 
It is thus readily seen that the effect of profile thick- 
ness is in many cases quite pronounced. The discrepancies 
are even more pronounced near the nose where the location 
of the forward stagnation point by the Birnbaum approxi- 
mation becomes very inaccurate, and infinite velocities 
occur at the mean line on flowing past the leading edge. 



5. Appendix: Structure of Uew Profile Forms 

In the ensuing development of new profile forms the 
computed source d is tr ibut i ons of the Kar man-Tr ef f t z pro- 
file are duly considered, since arbitrarily chosen source 
distributions would merely pr oduce " pr of i le formsthe ulti- 
mate shape of which it would be impossible to foresee. 
The only essential thing will he to so modify the form 
of a Karman— Tref f tz profile that the nose is given a 
fuller form, the location of the maximum profile thickness 
is changed, or a different mean— line form is prescribed; 
this means the application of the source distributions 
leading to the Joukowski profile (fig. 33a) and to the 
crescent (fig. 33c) with edge angle 6 and superposition 

by the simplest possible additive distributions f. 

s 

/ _ v N 

with a flow function Q, [ ~; ! . Then the stream func- 

\a • a « J 

tion of a symmetrical profile for the flow in x axis 
direction is given by: 




(98) 
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where 



tion 



z 



is a constant. The type of additive distribu— 
depends upon the desired profile shape. 



For a fuller nose, for example, a source and a sink of the 
same intensity is applied near the nose, or a doublet. 

As to the displacement of the point of maximum pro- 
file thickness, on Joukowski profiles, it is about one- 
fourth aft of the leading edge; on Karman-Tref f tz profiles 
it moves with increasing edge angle 5 closer toward the 
center. Suppose that a greater shift of maximum thick- 
ness toward profile center is desired. This displacement 



is favored "by the function f. 



interpreted as source 



superposition (fig. 10), which heretofore appeared solely 
as quadratic small term and which by itself has a crescent 
(fig. 33b) with angle 8=0 as streamline. The maximum 

thickness occurs at where the sources change to 



sinks. Selecting 



distribution follows at 



the total source 



Va' J " B.J 1 + I / a ' V 



where approximately |^for 



1 + 2 )+ i 5 JL 




(fig. 9). Putting 
file thickness 



Va »y 




a' J 

the point of maximum, pr o— 



(99) 



d 

Then the thickness is no longer computable from — and 

h a 
8, because — also exerts no effect upon it. There- 



fore, the thickness is prescribed at a selected point 
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D y 

'sVt'hat — "= 2 — . The s tr-eainl ine -r-epr e se nt ing a sym— 
a' a' 

metrical profile — curved profiles "being discounted ~ 
must therefore pass through the point with the coordi— 
*i ID 

nates — — and — — — , and carry the identification 
^, a 1 2 a 1 

S — s 0. Therewith equation (98) "becomes 

a' 7 X 

ID d n f X ! IDs 1 . / x , 1 D \ 

_ = _ t 4l , ; !+-h 5 --; 

2a' a V a » 2a'/ 2 Va ' 2a'/ 

h ,' Xl 1 D > 

+ — Q 4 V ( 10 °) 

a « \a» 2 a ' / 
D i. 

Besides the thickness — at — the point of maximum 
x m 



thickness — r and the edge angle 8 are given. The 
a * Xl 1 D N 

values of Q v ( — ; — — ) "being obtainable from figures 

V a ' 2 a « / d h 

15, 16, and 20, the two constants - and — themselves 

a a ' 

follow from equations (99) and (100) . The result, for 
cambered profiles, to which the known vortex distribu- 
tions are applied, is the flow function 

y y s d ' x y \ 1 _ / x y v 
a( 7^ a 1 a' a Va ' a ' y 2 5 Va ' a ' ^ 

+ <* 4 2 tan Pl | W ' i;l- > )- ¥ B (1|0)1 

a' Va'a'/ t- V a ' a ' ^ J 

+ 6 |^W 4 ^-L. -L^- W 4 (1;0) jj- (101) 

ys 

where, according to equation (53), the ordinate - — of 

a ' 

the plane of the singularities is given by 



4 
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- "(tan Pi j (1 ! (102) 

a 1 \ 2 a ' / V a' 2 / 



Since only linear term,? are involved, tan pj and \x 1 
no longer afford an accurate form of the median line of 
the profile (the form of the mean line remains for dis- 
appearing thickness only) as is readily seen on the first 
approximation in figure 21 (third profile). However, 
the camher of this median line can he established by pre— 

/ x y \ 

scribing a point P ( — ; — > 0 j as profile point and 



a ' a 



computing tan P 2 for this point from equation (101); 
whereas \i x is defined from the desired 0 angle of at- 
tack a 0 = —(Pi — ] . It is best to select the point 

V 4 J 

x y 
with the abscissa — =0 and the ordinate — = tan P 

1 d y f a ' D a ' 

+ — — - = — — + — , where — — can be given according to 

2 a 1 a 1 a ' a ' 

equation (100) or computed for a symmetrical profile by 
means of equation (98); tan P is the desired profile 
camber tan P ^ tan 

To illustrate: Take tan Pj = 0.08, m-i = 0.038 
and 5 = 16° from the NACA 4418 airfoil section and 
shift the position of the maximum thickness toward point 

x D x 

— = -0.2; while — = 0.32 remains for — =0. Then, 
a 1 a ' , a 1 

according to (99) and (lOO) , - *» 0.096 and — = 0.068. 

a 



Next is computedthe value of the flow function 



_i_x_ 

a ' V 



x 



from (101) and (102) for different abscissas and 

1. < 

{ * 



y 

several ordinates — — = 0.0, 0.05, 0.1, 0.15, by post— 

a 1 



om 



ing the values Q„ (— J — ' and V ( — ; — fr 

u \a ' a 1 / u \ a ' a ' y 

the figures 15 to 20 or tables 1 to ? and plotting the 
flow function against 

L y y y s 

a 1 a ' a 1 
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in figure 34. The values -—• for ■ ■ , • ■■ ■ = 0 are the 

■a' a' 7 X . 

$X 

points of the streamline — =0 of the looked— for 

a« V x 

profile. (The method of computing a profile from the 
singularities corresponds to that by Fuhrmann (reference 
6) and can he used for any arrangement of singularities.) 
For the example in question the profile of figure 35 is 
obtained. The nose is scarcely changed by the distribu- 
tion f it remains the same as on the Earman— 

Trefftz profile, whence its curvature circle is given 
with sufficient accuracy by equation (15) according to 
reference (3). The foremost profile point is by 
1 o / x N > 

- ^-from point A 2 1 ' — - = - 1.0 ) . The profile presents 
2a' V a 1 J 

a substantially different thickness curve from a normal 
Karman— Tref f t z profile and resembles those developed by 
Piercy, Piper, and Preston (reference 9) by conformal 
transf or mat i on . 

According to previous considerations, equation (59) 
can be retained for the flow along the y axis. The speeds 

— and — are lastly obtained by differentiation of 

(98), (101), and (59). The individual portions can be 
read off from figures 34 to 28. Since, in this instance, 
the profile form has been defined from the singularity 
distributions, the pressure distribution itself can be 
computed for these profiles also. Admittedly, a close 
approximation is obtained only for small a, since the 
streamline \Jr y is usually not exactly fulfilled. 



Translation by J. Vanier , 
National Advisory Committee 
for Aeronautics, 
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|_ o x/o 
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dx/a' 
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Quadra+ic 
term 
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d sr/a' 
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7 (Equafion 85). 
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Table 1. Flow function Q t of the source 

distribution f 1 (f/o')of the Joukowski Profile. 



o 




0 


0,05 


0,1 


0,15 


0,2 


0,3 




-0,75 


—0,191 


—0,210 


—0,221 


—0,231 


—0,242 


— 0,264 




-0,50 


—0,450 


—0,449 


—0,448 


—0,446 


—0,444 


— 0,417 


- 


-0,25 


—0,718 


—0,698 


—0,672 


—0,652 


—0,640 


— 0,605 




0 


— 1,028 


— 0,970 


—0,914 


—0,868 


—0,834 


— 0,764 




-0,25 


—1,220 


— 1,148 


— 1,083 


— 1',026 


—0,965 


— 0,860 




-0,50 


— 1,328 


— 1,223 


— 1,138 


— 1,061 


—0,991 


— 0,853 




-0,75 


— 1,162 


— 1,050 


— 0,964 


—0,892 


—0,824 


— 0,729 




-1,00 


0 


—0,318 


—0,376 


—0,403 


—0,414 


— 0,439 



Table 5. Flow function wA^j-,^) of the vortex 



distribution f i (l/a') of the S~ shaped mean line. 



\. V 














£\a' 
a' \^ 


0 


6,05 


0,1 


0,15 


0,2 


0,3 


0 


0 


0. 


0 


0 


0 


0 


±0,25 


T 0,127 


T 0,118 


T 0,108 


T 0,098 


=F 0,089 


T 0,068 


±0,50 


T 0,212 


T 0,192 


T 0,174 


T- 0,159 


T 0,145 


T 0,119 


±0,75 


T 0,238 


T- 0,216 


T 0,197 


T 0,179 


T 0,164 


T 0,139 


± 1,00 


T 0,165 


T 0,164 


=F 0,157 


T 0,148 


T 0,142 


=F 0,129 



Table 2. Flow -function Qn^' %)°* +hc source 
diE+nbu+ion f. ({/a') °f the-cresccnt 



\ V 

£ \o' 
«' \ 


0 


0,05 


* 0,1 ! 0,15 

1 


0,2 


o, - t . 


0 

±0,25 
± 0,50 
±0,75 
± 1,00 


—0,477 
—0,449 
—0,352 
—0,195 
0 


— 0,442 

— 0,405 

— 0,324 
—0,195 
—0,036 


— 0,420j — 0,401 
— 0,38Z*\— 0,369 
—0,304 —0,289 
—0,190—0,185 

— 0,055 ,—0,069 


— 0,385 
—0,350 

— 0,279 
—0,177 

— 0,080 


— 0,334 

— 0,312 

— 0,253 

— 0,177 

— 0,098 



Table 6. Flow-function Q 4 1-^; -^r) of the source 



distribution f 4 (l/a') ( F~igura 33). 



?-V 

a' X, 


0 


0,05 


0,1 


0,15 


0,2 


0,3 


0 


—0,339 


— 0,316 


—0,298 


—0,277 


—0,259 


— 0,231 


± 0,25 


—0,309 


—0,291 


— 0,274 


—0,256 


—0,239 


— 0,213 


±0,50 


—0,218 


—0,207 


—0,199 


—0,189 


— 0,180 


— 0,164 


± 0,75 


— 0,10 


— 0,100 


—0,102 


— 0,106 


—0,108 


— 0,109 


±1,0 


0 


— 0,018 


-0.030 


—0,042 


— 0,051 


— 0,060 



Table 3. Flow function Wt(%'< o') nf +he vortex 
distribution /„(!/«') of the flat plaTe. 



Table 7. Additive values to the Flow function 



X 

a 


\ -- 


0 


0,05 


0,1 


0,15 


0,2 


0,3 




- 0,75 


+ 0,75 


+ 0,78 


+0,80 


+0,82 


+0,84 




1-0,87 




-0,50 


+ 0,50 


+0,53 


+ 0,566 


+ 0,595 


+ 0,626 




-0,68 


1 


^0,25 


+ 0,25 


+0,280 


+ 0,319 


+ 0,364 


+ 0,407 




- 0,483 




0 


0 


+ 0,05 


+ 0,100 


+0,140 


+ 0,200 




-0,30 




-0,25 


— 0,25 


—0,195 


— 0,134 


—0,064 


0 




- 0,127 




-0,50 


— 0,50 


—0,419 


— 0,338 


—0,250 


—0,160 




-0,003 




-0,75 


— 0,75 


—0,628 


—0,509 


—0,385 


—0,270 




- 0,065 




-1,00 


— 1,00 


—0,570 


— 0,396 


—0,254 


—0,136 


+ 0,062 



(x y \ 
-J I of the vortex 

distribution / 3 (l/a') of the curved plate. 



\ V 
£ \a ; 
a' \^ 


0 


0,05 


0,1 


0,15 


0,2 


0,3 


0 


—0,250 


—0,203 


—0,158 


—0,114 




-0,070 




- 0,010 


±0,25 


—0,220 


—0,175 


—0,126 


— 0,089 




-0,048 




- 0,032 


±0,50 


—0,125 


—0,085 


—0,045 


—0,008 




1- 0,028 




1- 0,095 


±0,75 


+0,033 


+0,064 


+0,092 


+ 0,120 




-0,147 




- 0,196 


±1,00 


+ 0,25 


+0,255 


+ 0,265 


+ 0,275 




-0,288 




- 0,315 







X) 






X * 

V \ a' 
a' \ 


—0,95 


— 1,05 


— 1,10 


—0,95 


— 1,05 


— 1,10 


0 

0,025 
0,050 
0,10 


—0,582 
—0,559 
— 0,531 
—0,500 


0 

— 0,100 
—0,170 
—0,258 


0 

— 0,047 
—0,107 

— 0,204 


—0,009 

— 0,020 

— 0.028 

— ,042 


0 

—0,007 

—0,0145 

—0,0233 


0 

—0,004 
—0,012 
—0,022 










\ a a / 


X £ 


—0,95 


— 1,05 


— 1,10 


—0,95 


—1,05 


-1,10 


0 

0,025 
0,050 
0,10 


— 0,039 

— 0,040 

— 0,054 

— 0,080 


0 

—0,016 
—0,026 
—0,045 


0 

—0,009 
—0,019 
—0,037 


— 0,95 

— 0,82 

— 0,68 

— 0,47 


—0,424 
—0,374 
—0,329 
—0,241 


— 0,184 

— 0,170 

— 0,154 

— 0,102 



































W, 


la' ' a' i 






—0,95 


— 1,05 


1,10 


—0,95 


—1,05 


— 1,10 




-0,205 


+ 0.290 




h 0,334 




-0,195 




- 0,157 




- 0,139 




-0,209 


+ 0,290 




- 0,330 




- 0,182 




- 0,151 




- 0,138 




- 0,217 


+ 0,292 




-0,327 




- 0,176 




- 0,148 




-0,135 




-0,232 


+ 0,301 




- 0,327 




-0,171 




-0,142 




-0,131 
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Figs. 1,2,3,4,5,7, 8,9. 
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Figure 2. -Velocity jump at profile 
mean line. 
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Figure 1. -Circle transformed on a 
carved Joukowsky-prof ile 



yigure 3. -Division of speed on mean 

line in normal and 
tangential components. 
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Figure 5. -The additive formation 

K * a'/ 





Figure 4. -Mapping of a circle on a 

symmetrical Karman-Treff ts 
profile. 

Figure 8. -Simplification of integral yV 
of flow function. 
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Figure, 9. -Source distribution 'ig^re 7. -Identification of complex 
s A\ - flow potential in the 

f 5 (ar)of crescent plane of the section, 
for edge angle 6=1.0 
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Pigs. 6,20,29. 



Figure 6.- Circle mapped on S- shaped mean line. 
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Figure 20.- Flow function Q 4 (2 i a') of ' the additive source 
distribution f 4 (|) . 




Figure 29.- Gottinger airfoil section Ho. 624. 
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VorAt* distribution 1 f 3 (t/*)-yi-@of 



Figs. 10, 11, 12,13, 14, 15, 16. 





Figure 13— Identification of mean line 
— -i ' and axis AjAg of a Karman- 
; Treff ts profile. 
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Figure 14. -Solution of the four profile 
Figure lO.-Functions ffcCfe f 3 (4) and constants of a generalized 

^4 (4) ^ or superposed 
vortex distribution density and 
related forms of mean lines. 



Ka'rman-Treff ts profile. 
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Figure 15. -Flow function QiCjvfj) of 
» source distribution fj 

(4) of the Joukowski profile plotted 
«, against ordina^eX~for the different 
Figure ll.-Functions fi (feand f 5 (£) abscissas £ . * 

for the superposed source 
distribution density and related 
coat ours. 
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Figure 16. -Flow function Q5 (i^of 



I source distribution 

Figure 12. -Decomposition of y ordinate (if) of ^ crescent, 
in plane of section. 
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Figs. 17, 18, 19,21. 
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Figure 18. -Flow function W 3 (-§;£) of 



P 



Jiow function W 3 (■§! 
vortex distribution 



*H **» curved P la *« Plotted against 
Z. for different A- • 
a' B 



Figure 17. -Flow function Wg (-§;p of v/J-.--) 
| vortex distribution f 2 ,Q 

("a') of the flat plate plotted against 

7 for different x (VV 2 0.;^=+l,OO) . 
a W 



Figure 19. -Flow function W 4 (-fv^of 

I vortex distribution f 4 

("a 7 ) of the S-shaped mean plotted 

against JL for different X . > 
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Pigs. 22, 23,33,34, 35. 
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Figure 22. -Example; 

approximation of 
three NACA airfoil 
sections. 




Figure 23. -Speed in vicinity of a 
profile. 



Upper surface 

d Profiles 





Figure 33. -Profile forms obtained 

by the individual source 
distributions in parallel flow along 
x axis. 
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Figure 35. -Hew profile forms by 

pre scribed special 
singularity distributions. 



Figure 34. -Example; the valuf of flow 
function %£h on individual 
points-^- for different ordinatea-^ by 
given singularity distributions. 



Figure 25. -Orthogonal system of velocities induced by Figure 24. -Orthogonal system of velocities induced by 
the source distribution fg (JL) of the crescent. the source distribution fj (|,) of the 

a ' Jowkowsky profile. 




Figure 30. -Speeds in field of 

Gottinger airfoil 
section No. 624; section -§-=-0,5 
parallel to y axis.Qmmission of 
profile thickness affords the 
dashed curve .disregard of 
quadrative term dSx / the dash- 
dot curve. dx ' & 
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Figure 31. -Speeds in field of 

Gottinger airfoil 
section No. 624; section -f?=+ 0,5 
parallel to y axis. 
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Figure 32. -Speeds in field of 

Gottinger airfoil 

section No. 624; sec ti on -l t 5 J* 

parallel to y axis. ~ o<j 
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